Noise spectrum of quantum transport through quantum dots: 
a combined effect of non-Markovian and cotunneUng processes 
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Based on our recently developed quantum transport theory in term of an exact master equation, 
the corresponding particle-number resolved (n-resolved) master equation and the related shot noise 
spectrum formahsm covering the fuU frequency range are constructed. We demonstrate that the 
noise spectra of transport current through single quantum dot and double quantum dots show 
characteristic steps and/or peak-dips in different tunneling frequency regimes through tuning the 
applied bias voltage and/or gate voltage at low temperatures. The peak-dips crossing the tunneling 
resonant frequencies is a combination effect of non-Markovian and cotunneling processes. These 
voltage-dependent tunneling resonance characteristics can be utilized to effectively modulate the 
internal Rabi resonance signature in the noise spectrum. 

PACS numbers: 72.10.Bg; 05.40.-a 



I. INTRODUCTION 

Shot noise of non-equilibrium current fluctuations con- 
tains rich information beyond the average current 
It has stimulated great interest in recent years both 
in theory and experiment .'^"l"^ The evaluation of shot 
noise depends on the development of non-equilibrium 
quantum transport theory. Conventional approaches to 
quantum transport include the Landauer-Biittiker scat- 
tering matrix theoryj^^ the non-equilibrium Green's 
function formalism /^i^^ and the real-time diagrammatic 
techniqueJ^'-iS Many interesting problems such as super- 
and sub-Poissonian noises and their origins have been 
addressed i^iiSiiSS However, most of the investigations fo- 
cused only on the zero- or low- frequency regions. The 
evaluation of noise spectrum over the full frequency range 
is largely lacking and remains challenging with these ap- 
proaches. Another commonly used method in quantum 
transport is the quantum master equation approach. In a 
certain sense it is simpler and more straightforward than 
the approaches mentioned above.— It has the advantage 
of generality for different scattering processes being han- 
dled in a unified manner, even the transient dynamics 
under the time-dependent bias voltage<^"— 

Quantum transport theory based on master equation 
has been developed rapidly and studied extensively in re- 
cent years However, most of work were based 
on a certain perturbative master equation description, 
which is valid mainly in Markovian dynamics dominated 
regime with sequential tunneling processes involving only 
one electron tunneling events at a timeji"— '^'^i^^i^^ Very 
recently, on the basis of Feynman- Vernon influence func- 
tional technique, one of us has obtained an exact mas- 
ter equation for nanostructures,— from which we have 
also established an exact non-equilibrium quantum trans- 
port theory for noninteracting electronic systemsi^ This 



new non-equilibrium quantum transport theory is appli- 
cable for arbitrary voltage, arbitrary temperature and 
arbitrary system-reservoir coupling strength. Quantum 
transport based on this exact treatment can deal with 
all the non-Markovian tunneling processes accompanied 
with not only sequential tunneling but also cotunneling 
involving two or more electron tunneling events at the 
same timei ^^i^° 

In this paper, we will construct the corresponding mas- 
ter equation for the reduced density matrix conditioned 
on the electron number passed through the tunnel junc- 
tion, i.e. the particle-number resolved (n-resolved) mas- 
ter equation. Together with the MacDonald's formula, 
we give the description for non-Markovian noise spec- 
trum over the full frequency range. All relevant formu- 
lations will be summarized and developed in Sec. |ll] to- 
gether with the detailed numerical method presented in 
Appendix. This formalism is certainly applicable to ar- 
bitrary bias voltage, temperature and system-reservoir 
coupling strength for all the non-Markovian processes. 
In Sec, mil based on the present formulation, we study 
the quantum transport setup for the shot noise of current 
through two model systems. The first model is a single 
one-level quantum dot system, operated in the conven- 
tional sequential tunneling regime when only the average 
current is concerned. From the noise spectrum, however, 
we find that there emerge peak-dips in the noise spec- 
trum as the temperature decreases. Another example is 
a coupled double quantum dots system, with one level in 
each dot. The resulting noise spectrum shows both the 
coherent Rabi resonance of the dots system and the ex- 
ternal field-dependent tunneling characteristic resonant 
structure. We observe that the coherent Rabi resonance 
is insensitive to the temperature, while the external field- 
dependent tunneling resonant structure has the very in- 
teresting temperature-dependent phenomena. In partic- 
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ular, the external field-dependent tunneling resonance 
characteristics is manifested with either a peak-dip or a 
step at low temperature, depending on the corresponding 
tunneling regime defined by the applied bias and / or gate 
voltage. Consequently, by varying the applied voltage 
one can tune the tunneling resonance on top of the Rabi 
resonance, thus effectively modulate the coherent signa- 
ture in noise spectrum. Finally, a summary is presented 
in Sec.HVl 



II. 



n-RESOLVED MASTER EQUATION AND 
NOISE SPECTRUM 



A. Exact master equation and transient transport 
current 

Let us start with a brief outline of the recently de- 
veloped exact master equation and the non-equilibrium 
quantum transport theoryj^i^ The quantum transport 
studied with a nano-device consists of such as quantum 
dots interacting with the electron reservoirs (source and 
drain). The total Hamiltonian of the system has three 
parts: Ht = Hs + Hb + H' . The central system of 
quantum dots part assumes Hs = ^^^iJ.i^O'liO'i^, where 
c-jj, i^fj.) is the creation (annihilation) operator of electron 
at the specified state of energy e^^ = e^, while the off- 
diagonal 6^1/ is the hopping integral. The Hamiltonian of 
free-electron reservoirs reads Hb — J2ak ^akcl^f^Cak, with 
cjjj, being the creation operator for the electron in the 
reservoir a = L,R. The system- reservoirs coupling for 
electron tunneling is given by 



H' 



H.c. 



(1) 



akfi 



where Vakf^ is the tunneling coefficient or the tunnel- 
ing amplitude between the system and the reservoir a. 
The electron-electron interaction has been ignored here. 
Then the reduced density matrix of the central system 
p{t) = tTB[pT{t)], i.e., the trace of total density matrix 
over reservoirs degrees of freedom, is governed by the ex- 
act master equatio: 
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p(t) = -i[Hs{t),p{t)] -^-fc,^,^{t)[al,a^p{t) + p(t)a^ 
- "^lat^y {t) (^ala^,p{t) ip(t)a|,a^ - a,,p{t)c 



(2) 



where the first term represents the Liouville equation (the 
unitary part) of the system with the renormalized Hamil- 
tonian Hs{t) — 'Yl^iu^t^vCL^^cLi' and the renormalized en- 
ergy levels e{t) = e{t) -^Ec ['*«(0 - H^l(<)] • The other 
two terms (the non-unitary part) in Eq. ^ describe the 
dissipation and decoherence dynamics of the system with 



the time-dependent dissipation and fluctuation coeffi- 
cients 7a (i) = ^[Ka{t)+Kl{t)] and 7a (t) = Xa{t)+Xl^{t), 
where Ka{t) and Xa{t) are given by, 

Ka{t)=l dTgait,T)u{T)[u{t)]''^ , (3a) 

Aa(t) = / dT{gait,T)v{T) - ga{t,T)u{T)} ~ Kait)v(t). 
•I to 

(3b) 

and u{t) and v{t) are indeed directly related to the 
Keldysh's non-equilibrium Green functions and they 
obey the following dissipation-fluctuation integrodiffer- 
ential (Dyson) equations:— 

u{T)+ie{T)u{T)+Y, f dT'g^{T,T')u{T')^0, (4a) 

a •'to 

v{t) + ie{T)v{T) + ^ / dT'ga(T, t')v{t') 

a ■'to 

= f dT'ga{T,T')u{r'), (4b) 

a ''to 

subject to the boundary conditions u{to) — 1, and 
v{to) = 0. The integral kernels in the above 
equations, ga{T,T') and ^^(t, r'), are defined by 

andga^^(r,r') = 



ku^ 



T,k^akf,V*^^fa{eak)e "'""'^^ , which depict all the 
non-Markovian memory processes of electrons through 
the tunnelings between the dot system and the leads, 
and fa{eak) = l/(e'3°(^°fe-^°) - 1) is the Fermi distri- 
bution function of the lead a in the initial equilibrium 
state. Introducing the spectral density functions of the a- 
lead coupled with the dot system: Taiuj) = {ra^^(w) = 
'^^^k ^aktxV*f^,^S{u! — iak)}, the memory kernels are re- 
duced to 



/duj 
— /a(a;)ra(a;)e 



(5a) 

tu(t~T)^ (5b) 



Through out this work, we have set e = h = 1. 

The transient current of the a-lead is obtained as2^ 

/a(i) = -2Re [ dTi:Ag„{t - t)v{t) - - ^Mt) 

•I to ^ 

+ ga(i,r)M(T)e(to)w^(t)}, (6) 

with u{t) = u^{t — t) and Q{t) being the single particle 
reduced density matrix, Qfiy{t) = tr^ [aj^ a^p(t)] , satis- 
fying Q{t) ~ v{t) + u{t) Q{to)u^t) . The above current 
expression can easily recover the Landuer-Biittiker for- 
mula obtained by scattering approach and the Keldysh's 
non-equilibrium Green function formulation, as detailed 
in Ref.lM 
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B. n- resolved quantum master equation 

Before constructing the n-resolved master equation, 
we should point out that the master equation of Eq. ([2]) 
was derived exactly for noninteracting quantum dot sys- 
tems so that all the electron tunneling processes in such 
systems, including sequential tunneling and cotunneling 
processes together with all the non-Markovian memory 
effects, are fully taken into account. For noninteracting 
systems, the cotunneling processes refer to the tunnel- 
ings involving two and more electrons through the junc- 
tions simultaneously. In the conventional master equa- 
tion in terms of perturbation expansion, one may expect 
that the sequential tunneling and cotunneling processes 
should be described separately by one pair and two or 
more pairs of electron creation and annihilation opera- 
tors, respectively, in the non-unitary part of the mas- 
ter equation. However, for noninteracting systems, the 
master equation derived from the perturbation expan- 
sion up to the second order of the tunneling coefficients, 
i.e. Vak^i in Eq. ([1}, has indeed the same operator struc- 
ture as that of the exact master equation of Eq. ([2]).^^ 
The only difference between the second-order (2"'^-order) 
perturbative master equation and the exact master equa- 
tion is manifested in the determination of the dissipation 
and fluctuation coefficients in the corresponding master 
equations. In the exact master equation of Eq. these 
dissipation and fluctuation coefficients, i.e., "fafivit) and 
"fafj,u(t), are determined by the Green functions u{t) and 
v{t) which satisfy the exact dissipation-fluctuation inte- 
grodifferential equations of Eq. The non-Markovian 
memory effect and the cotunneling processes described 
in the exact master equation are fully characterized by 
the non-local integral kernels in Eq. (j?]), i.e. the Dyson 
equations through the iteration to all orders of the tun- 
neling coefficients. Truncating Eq. Q to the 2"'^-order 
perturbation of the tunneling coefficients, all these time- 
dependent dissipation and fluctuation coefficients are re- 
duced to the coefficients in the 2"''-order perturbative 



master equation that can be derived directly from the 
quantum Lioulille equation in the perturbation expan- 
sion approach. Since the 2"'*-order perturbative master 
equation can only describe the sequential tunneling and 
Markovian process, it indicates that the cotunneling pro- 
cesses and non-Markovian effect in noninteracting sys- 
tems are fully determined by the time-dependent dissi- 
pation and fluctuation coefficients from the solution of 
Eq. ^ , rather than the operator structure of the master 
equation. The operator structure of the 2'^'^-order pertur- 
bative master equation is exactly the same as that of the 
exact master equation for noninteracting systems. The 
detailed derivation of this connection between the exact 
master equation and the 2"'^-order perturbative master 
equation has been given in Ref.[28l 

Based on such a connection between the exact mas- 
ter equation and the 2"'^-order perturbative master equa- 
tion we just discussed above, we can construct the cor- 
responding n-resolved master equation from Eq. ([2]) for 
the conditioned reduced density operator from the 

n-resolved 2"'^-order perturbative master equation that 
has been derived explicitly from the 2"'^-order perturba- 
tive master equation. ^("q) ig conditioned on the reg- 
istered number Ua of electrons having passed through 
the tunnel junction between the specified a-lead and the 
central dots system. As it has been shown in Ref.[2^, the 
resulting n-resolved 2"'^-order perturbative master equa- 
tion only modifies the corresponding non-unitary opera- 
tor structure of the 2"'^-order perturbative master equa- 
tion, namely, it replaces only the registered jump terms 
aj^pa^, and a^paj^, associating with the specified a-lead 
in the master equation by ajjp(""+^^ay and a^p(""^^)aj,, 

respectively. Other terms remain the same but for 
Since the exact master equation for noninteracting sys- 
tems has the same operator structure of the 2"'^-order 
perturbative master equation, by analogy, we obtain the 
corresponding n-resolved master equation of the exact 
master equation of Eq. 



\t)a,. (7) 



J2 [7.M^(i) +27aM.W]«^/"°''^W«|. -T.^c.,AtKp^""^'^i 



for keeping track the number of electrons tunneling only and 
through one specific lead, e.g., the a-lead with a' ^ a. 



4 



for conditionally keeping track the number of electrons 
tunneling through both the right and the left leads. Here 
we omit the time index in the time-dependent coefficients 
and have also defined jf^i, = J2a=L B.'^o'f^'^- Both the n- 
resolved master equations of Eq. ([7]) and Eq. ([SJ will be 
used to derive the noise spectrum of an individual lead 
and the cross-correlation noise spectrum. 

The same as the connection between the exact mas- 
ter equation and the 2"'^-order perturbative master equa- 
tion discussed in the beginning of this subsection, the n- 
resolved master equation, Eq. ([7]) and Eq. ([5]), also has 
the same operator structure as that of the n-resolved 
2"'^-order perturbative master equation. •^^ The only dif- 
ference is the dissipation and fluctuation coefficients in 
these master equations. The dissipation and fluctua- 
tion coefficients in the rt-resolved 2"'^-order perturba- 
tive master equation involve only sequential tunneling 
process and the Markovian dynamicsi^ However, the 
time-dependent dissipation and fluctuation coefficients in 
Eq. ([7]) and Eq. ([5]) that are determined by the exact non- 
equilibirum Green functions of Eq. Q for noninteract- 
ing systems have taken into account simultaneously the 
sequential tunneling and cotunneling processes, so does 
for the non-Markovian dynamics. One may ask why the 
states such as ^(„^_i,„^+i) ^(„^+i,„^_i)^ 

etc. which are expected to be directly related to the co- 
tunneling processes for the electron registered detector 
in the perturtive expansion do not occur in Eq. ([7]) and 
Eq. ([8]). This must be the same reason as in the case of 
the master equation itself. In terms of the perturbatve 
expansion, the master equation up to the higher order 
contains the dissipation and fluctuation terms involving 
two and more pairs of electron operators sandwiching the 
reduced density matrix that describe the cotunneling pro- 
cesses. But the exact master equation does not contain 
such terms. This is because the cotunneling processes 
has been all switched into the dissipation and fluctua- 
tion coefficients determined by the exact Dyson equa- 
tions. Although we have not provided a rigorous deriva- 
tion of Eq. ([7]) and Eq. it should be the same reason 
that the n-resolved master equation, Eq. ([7]) and Eq. ([8]), 
can maintain the simple operator form because the time- 
dependent dissipation and fluctuation coefficients deter- 
mined by the exact nonequilibrium Green functions can 



I 

address all the sequential tunneling and cotunneling pro- 
cesses as the non-Markovian memory dynamics. 

As a justification and also a self-consistent check, 
from the n-resolved master equation of Eq. ([7]) with the 
identity p{t) = ^„ p*^"°^(i), we can easily reproduce 
the exact master equation of Eq. ^ . The exact cur- 
rent expression of Eq. ([6]) can also be reproduced from 
the n-resolved master equation of Eq. ([7]) as follows. 
With the knowledge of p'^"°^(t), the tunneling electrons 
distribution can be readily evaluated via P{na,t) = 
Xr5p("°)(t). This is the key quantity for full counting 
statisticsj^i^'^ The m*^-moment of transport is just 
(n™(i)) = from which all transport 

properties can be obtained. For instance, the measured 
current which is related to the rate of the first moment, 
is given by Ia{t) = —■^{na{t)). Using the n-resolved 
master equation of Eq. (|7]) , it is also easy to reproduce 
the exact expression of the transient transport current 
of Eq. ([5]) . The capability of recovering the exact master 
equation of Eq. ^ and the exact transport current for- 
mula of Eq. © from Eq. ([7]) ensures that the n-resolved 
master equations we constructed here is the current one 
for noninteracting systems. Thus the current-current cor- 
relation noise spectrum which is related to the 2"'^-order 
cumulant can be directly evaluated now from Eq. ([7]) and 
Eq. ([8|), which will be presented in the following. 



C. Noise spectrum expression 

We consider now the current noise spectrum, S{uj) ~ 
,!^{{5I{t)5I{0))s}, i.e., the full Fourier transformation 
(denoted by ^) of the fluctuating current-current corre- 
lation function that is symmetrized. For the total circuit 
current I(t) = alL(t) — bIfi{t), which is typically the mea- 
sured quantity in most experiments^ with the coefficients 
satisfying a + b = 1 related to the symmetry of the trans- 
port setup (e.g., junction capacitances), the circuit noise 
spectrum is S{uj) = a'^SL{u}) + b'^SR{uj) -2abSLRiuj)-^^ 

The noise spectrum at individual lead a = L or i? is 
given by the MacDonald's formula, Sai^j) = S'qq(w) = 
2uj J^dt sm{ujt)§[{nl{t)) - {Lt)'^]. With the help of 
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Eq. (O, the involved quantity {n\{t)) satisfies 

j^{nl{t)) = 2tr [{%+jc)JVc+%Afa] 

+ tr [(7„ +7„)£.-7„^] , (9) 

where ^ = 1 — ^ is the single hole particle reduced den- 
sity matrix, while Afa^^i^ = trs[aJ,a^iV"] and Afa^j,!^ = 
tTsia.fiCilN"'] are related to the number operator N°'{t) = 
X]n„ ^aP^^^Ht)- It satisfies 

= MN"^ + tr [(7„ -t- 7„)e] + tr [7^5] , (10) 

with M. being the generator of Eq. ^ that is recast as 
p = Mp. The initial conditions to Eqs. ([9]) and (|T0)) are 
the stationary states without the electron tunneling, i.e., 
iV„(0) ^ and p(0) = p^'K We have 

Sc^iw) = Awlm^ {tr[(7a + 7a)A/'a + laJ^a] } 

-H2Reif{tr[(7,+7„)e''-7„r*]}, (H) 

where ^{f{t)} = dte''^*f{t) denotes the Laplace 
transformation. 

The cross correlation noise spectrum is given byi^i^i^i 
Slr(.uj) = i^{(«L(t)««(0))s + {6lRit)SlLms} = 
2L,j;^dtsm{ijt)^[{NL{t)Nnit)) - (It)^], where 
{NL{t)NR{t)) = trE„,„«"iW™HO- Again, 
with the help of Eq. ([8]), we obtain 

Slb.{lo) = 2cjlmif |tr[(7i + jL)AfB. + IlMr] 

+ tr [(7fl + 7Ji) ATi + 7i?.A^i] } ■ (12) 

We thus have completed the expressions of the current 
noise spectrum, i.e., Eq. (fTTj) and Eq. ([T2|). which are the 
main results of the present work. They can be further 
written in the compacted formalism as 

5'aa'(w) = 2wlm^|tr[(7a +Ja)J^a' + laJ^a'] 

+ tr[(7Q' + 7a')A/'a + 7a'A/'a] | 

+ 2J„„,Re^{tr[(7, + 7„)e'' - 70^*]}- 

(13) 

It describes both the sequential tunneling and cotun- 
neling processes together with the non-Markovian mem- 
ory effect involved in the time-dependent dissipation and 
fluctuation coefficients Jait) and 7a (t) that also deter- 
mine the master equation of Eq. ^ and the n-resolved 
master equation of Eq. ([7]) . Wc will evaluate explicitly the 
noise spectrum using Eq. in the coming section with 
two widely adopted quantum transport model systems. 
Since the time-dependence dissipation and fluctuation co- 
efficients are determined by the functions u{t) and v{t) 
which are governed by Eq. Q , the key to understand 



the cotunneling processes and the non-Markovian effect 
in the transient transport current and the noise spec- 
trum is to solve exactly the integrodifferential equations 
of Eq. The detailed numerical method is presented 
in Appendix based on a parameterization scheme. 

As it is shown in Ref.Hl, u{t) and v{t) accounts for all 
orders in perturbative expansion. If one solves Eq. (|4]) up 
to the second-order of the system-reservoir couplings, the 
resulting master equation recovers the 2"'^-order master 
equation in the perturbation theory that has been stud- 
ied widely in the literature, see for examples Refs. 6, 7,il. 
The usual Markovian approximation is obtained under 
both the wide band limit Ta{Lu) = Fq, with high tem- 
perature, which leads to gQ(^ — r) — > ^S{t — r) and 
ga{t — r) — >■ ^faS{t — r). Then all the time-dependent 
coefficients in the master equation become constants, 
e.g., 7ct(i) = and 7a(i) = —fcJ'a, which recovers the 
n-resolved Markovian master equation in Ref.lsH and also 
reproduces the solution given in Ref.|3llfor resonant-level 
model. 

Usually one believes that in the wideband limit and the 
high bias voltage regime, the non-Markovian effect van- 
ishes and the Markovian master equation works. How- 
ever, the non-Markovian memory effect is very compli- 
cated, it relates to several factors together, such as the 
temperature, the bias voltage, the spectral bandwidth, 
and the system-reservoir coupling ^rength as details in 
Ref.ll^ and in Appendix. In Ref.,23!, we showed that in 
the WBL, the large bias voltage indeed weakens the non- 
Markovian effect, but does not lead to a complete Marko- 
vian limit. At this regime, the non-Markovian effect, 
which arise completely from the low temperature and/or 
strong system-reservoir coupling strength, may not be 
clearly manifested in the transport current itself and the 
reduced density matrix, ^'^ but it becomes significant in 
the noise spectrum of the current-current fluctuation^^ii^ 
as shown in the next section. 



III. NUMERICAL DEMONSTRATIONS 

A. Single quantum dot system 

For simplicity, we consider in this section the situ- 
ation of an energy-independent spectral density of the 
leads, namely, the flat band of ra(a;) = Fq. and denote 
F Tl + Tji. We set the symmetric junction capac- 
itances as a = 6 = 0.5 and also the symmetric bias 
oi PL = —PR = eV/2 in the following studies. Let 
us start with a single one-level dot model system, with 
Hs = ea^a. The average current of this system can be 
evaluated as / = TlTr J (flL^l+^g/z)^ ■ Despite of its 
simplicity, such a single quantum dot system has many 
interesting physical phenomena and has been widely used 
as a single electron transistor. Here we operate the sys- 
tem in the large voltage regime, p^ > e > pn, which is 
also considered as the conventional sequential tunneling 
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FIG. 1: Noise spectrum (in unit of e^V/h) for current 
transport through single quantum dot, where e = F, with 
Tl = Fij = 0.5 F, under the bias oi V = lOF at the speci- 
fied temperatures. The peak-dips crossing the tunnehng res- 
onance frequencies, ujc, = \e — fia\ with a — L and R, are 
indicated by arrows. The inserting figures in (a) and (b) give 
a comparison between the exact solution with the 2"'*-order 
perturbation result at the temperature ksT = O.IF. 



regime if the average current is concerned. However, we 
will see below that in this case a striking non-Markovian 
character, manifested by a peak-dip, will emerge in the 
current noise spectrum. 

Figure [U shows the current noise spectrums for 
auto-correlations [Fig.[TJa) and (b)], cross-correlation 
[Fig.[TJc)], and circuit current correlation [Fig.[TJd)], at 
different temperatures. In the low or high frequency re- 
gion, the noise spectrum is nearly independent of tem- 
perature, despite of the fact that the stationary cur- 
rent increases as temperature decreases. In the large 
frequency regime, the noise spectrum at the individual 
lead Sa{ci>) r^, while that of cross-lead S'Lfl(i^) — > 
0. More interestingly we observe that a peak-dip fea- 
ture is developed in the noise spectrum of current auto- 
and cross-correlations, crossing each resonant frequency 
Wq = |e — /IqI (indicated by an arrow in Fig.[T]) as the 
temperature decreasing. 

The above observations can be understood as follows. 
As the zero-frequency noise and also the stationary cur- 
rent are concerned, the transport system studied here 
with /ii > £ > is dominated by sequential tunneling 
processes. In the high frequency (w ^ |e — fia\) region, 
we have Slr{i^) ~^ 0, due to the electron correlation be- 
tween different leads is yet to be established. However, 
Sa{tL)) — Fq coming from the fact that the fluctuations 
arises mainly from the reservoir background. Therefore, 
we cannot find the system-associated structure in the 
noise spectrum at both limits. However, the shot noise 
over the full frequency range describes various tunnel- 
ings associated with energy emissions and absorbtions at 
different detection frequency u. Therefore it must man- 
ifest the energy structure dependence of the dot system 
through the applied bias voltage and temperature. As 



one can see from Fig.[T] the noise spectrum at relative 
high temperature shows a smooth energy-structure de- 
pendence near the resonant frequency Wq = |e— /icj. This 
result is consistent with that of Refs.|3 and[l3l. in which 
one used the second-order time-nonlocal master equa- 
tion which contains very little memory when fc^T > T 
and describes only the sequential tunneling process where 
the Markovian dynamics is dominated. However, non- 
Markovian dynamics should play an important role in the 
noise spectrum at low temperature {kgT <C F) where co- 
tunneling processes should become significant, especially 
around the resonant frequency where the electron can 
dramatically tunnel forth and back between the leads 
and the dot near the Fermi surface of the leads. The ob- 
served peak-dip crossing the resonant frequency in S'a(w) 
[Fig.mja) and (b)] and S'Lfl(i^) [Fig.fTJc)] is indeed such 
a non-Markovain result as a combination of sequential 
tunneling and cotunneling processes at low temperature. 
In particular, the positive peak-dip feature in the cross- 
correlation fluctuations spectrum must come from the 
cross-lead cotunneling contribution. 

The above peak-dips crossing the resonant frequencies 
in the noise spectrum have not been shown in the pre- 
vious studies-ii^-^i^ii^iSirlg ^his is because in Refs.i-H 
the Markovian treatment has been used, while Refs.|j 
and [13 considered only the sequential tunneling process 
at relative high temperature with the 2"'^-order mas- 
ter equation approach which is also Markovian dynam- 
ics dominated. Refs.fl2| used Biittiker's scattering ma- 
trix approach which is exact at zero frequency but con- 
tains approximation that covers mainly the Markovian 
dynamics at finite frequencyF^^i^ To make a compari- 
son, we calculate the noise spectrum with 2'^'^-order per- 
turbation approximation at low temperature. The re- 
sult is inserted in Fig.[IJa) and (b). It shows that the 
2"'^-order perturbation noise spectrum shows a dip only 
but not a peak-dip. As it is well-known the 2"'^-order 
perturbation approximation is invalid at low tempera- 
ture. Very recently, using the nonperturbative hierar- 
chical equation of motion treatment^ this peak-dip fea- 
ture in the noise spectrum is also observed. Fig.[T] also 
clearly shows that with the temperature increasing, the 
aforementioned non-Markovian effect, i.e. the peak-dips 
in the noise spectrum, is diminished. We thus conclude 
that the observed peak-dip feature in the noise spectrum 
is a non-Markovian memory effect through various tun- 
neling processes, including both the sequential tunneling 
and cotunneling events. 



B. Double quantum dots 

Consider now the transient transport through a sys- 
tem of two coupled quantum dots, described by Hs = 
J2fi=i r ^M^i'^M + ^{djdr + d\.di). This system has been 
studied widely as a charge qubit^'i^i^'^S and has also 
been proposed as an alternative detector of a charge 
qubit.'^^ The intrinsic coherent Rabi frequency A is the 
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FIG. 2: Shot noise (in unit of e^F/ft) of a double-dots system, 
where ei = Er ~ and = 5F, with Tl = F_r = 0.5 F at 
the specified temperatures. Left panels: The noise spectrum 
as a function of frequency, under the bias voltage eV = 20 F. 
The tunneling resonances of the peak-dip chacteristics occurs 
at ij± — \e± — l-ial ~ 5F and 15 F and the Rabi resonance 
at A = 10 F, indicated by the arrows individually. Right 
panels: The shot noise at the Rabi frequency as function of 
bias voltage. The double-resonance uj — A = \e± — fJ.a\ occurs 
at both eV = lOF and eV = 30 F. The former satisfies also 
the resonant transport conditions fiL = £+ and hr = e_, but 
the latter is of /iL > e± > 



energy difference between eigenstates (e±), e.g., A = 
e+ — e_ = 2n for the degenerate double-dots system 
considered for demonstrations. It is known^i^S that the 
Rabi coherence of the central system shows a dip at 
a; = A in the auto-correlation noise spectrum £'0(0;), 
as can be seen in Fig.[2fa). But, it appears a peak in 
the cross-correlation noise spectrum 5*^^(0;), as shown 
in Fig. [21b), while remains a dip in the total circuit noise 
spectrum, for symmetric junction capacitances, as de- 
picted in Fig.[2Ic). The above Rabi coherence signa- 
tures are nearly independent of the temperature. Phys- 
ically, the Rabi coherence is intrinsic. Therefore, it 
can be extracted even in the weak system-reservoir cou- 
pling regime, where the 2"'^-order master equation is 
applicable Besides the coherent signals of Rabi fre- 
quency in Fig.[5][a)-(c) where fj,L > £± > ^J■R, the ex- 
pected peak-dips of non-Markovian characteristics occur 
at tunneling resonance of uj± = \e± — ^ia\ at low temper- 
ature. As temperature increases, this peak-dip feature is 
diminished, just the same as that in the single quantum 
dot case described earlier. 

Figure [2I^d)-(f) depict the noise spectrum at Rabi res- 
onance u! = A versus the applied bias voltage V that 
monitors the tunneling resonance of cj = \e± — /icl. 
The double-resonance of w — A — \e± — fj,a\ occurs 
at y = lOr and V — 30r for the present system in 
study. Apparently, the bias voltage dependence of the 
noise S{A) at Rabi frequency is changed dramatically 
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FIG. 3: Frequency-dependent noise spectrum of the left lead 
for the double-dots transport system with different tunneling 
regimes through gate voltage Vg, i.e., ei — Sr ~ Va, at low 
temperature ksT = 0.1 F with the applied bias voltage eV = 
5F. The other parameters are the same as in Fig.[2l 



at low temperature, especially at tunneling resonance, 
as an expected non-Markovian effect. In particular, at 
V = 30 r, both the two dot levels lie within the applied 
bias window, i.e., > e± > fiR. The resulting peak- 
dip tunneling resonance shows up clearly in Fig.[2l^d)-(f). 
The behavior of S{uj — A) around V = 10 F is an inter- 
play between the Rabi resonance and the lead-dot tun- 
neling resonance. 

The frequency-dependent noise spectra demonstrated 
in Fig.[T] and Fig.[2][a)-(c) are operated in the so-called 
large bias region of fiL > £, e± > Mfli where the tunneling 
resonance shows the peak-dip feature at low temperature, 
as a combination effect of the sequential tunneling and 
cotunneling events. Since the Rabi resonance spectral 
profile does not depend on the applied gate voltage for 
its overall intensity, we shall also examine the tunneling 
resonance behavior for different set up of the applied gate 
voltage to the dot energy levels. Without loss the gener- 
ality, we take a fixed bias voltage, say V = 5T, and apply 
several different gate voltages to the dot levels to achieve 
different tunneling scenarios. The results are summa- 
rized in Fig. [3] Fig.[31Ja) and (c) describe the case of both 
the two dot levels lying outside the bias voltage window. 
The corresponding tunneling resonance shows only with 
a step shape. In contrast, if one of the dot levels, i.e., 
the level £_ and e+ in Fig.[31[b) and (d), respectively, lies 
inside the the bias voltage window, the corresponding 
tunneling resonance spectrum generates a peak-dip. As 
it has been understand, for the dot level lying outside the 
bias voltage window (i.e. either e > fJ-L,R or e < hl,r), 
the sequential tunneling through this level is largely sup- 
pressed and the transport is dominated by cotunneling 
events. So the cotunneling transport alone only shows a 
step near the resonance frequency. The peak-dip occurs 
only for > e > fj,R at low temperature, where both 
the sequential tunneling and cotunneling participate in 
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the non-Markovian dynamics. 
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IV. SUMMARY 

In summary, based on the exact master equation 
for nonequilibrium transport, we constructed the corre- 
sponding particle-number resolved master equation, from 
which we have also established the formalism for the 
noise spectrum in the full frequency range. This new 
formulism is applicable to arbitrary bias voltage, tem- 
perature and system-reservoir coupling strength for all 
the non-Markovian processes in various nanoelectronic 
systems where the electron-electron interaction has been 
ignored. We applied this formalism to two widely studied 
transport model systems, i.e, electron transport through 
single quantum dot with a resonant level and double cou- 
pled quantum dots containing one energy level in each 
dot contacted the electron reservoirs, respectively. We 
showed different tunneling characteristics in the noise 
spectrum in both the large and small bias voltage regime 
at different temperature. In particular, we found peak- 
dips in the noise spectrum crossing the tunneling reso- 
nant frequencies which are defined individually by the 
energy difference between the applied chemical poten- 
tials and the dots energy levels. This peak-dip feature 
has not been observed in previous theoretical works. The 
characteristic peak-dip in the current noise spectra is a 
non-Markovian effect at low temperature involving both 
sequential tunneling and cotunneling. As the tempera- 
ture increasing ((fc^T > F), peak-dip profile in the noise 
spectrum is diminished. In contrast with the aforemen- 
tioned external bias voltage regulated resonant character- 
istics, the internal coherent Rabi oscillation signal in the 
double dot system is rather independent of temperature. 
The coherent Rabi oscillation results in just a normal dip 
profile in the auto-correlation noise spectrum but a peak 
in the cross-correlation spectrum. We expect that these 
characteristics in the current noise spectra we found here 
would be tested in experiments or in other theoretical 
calculations. 

We should also point out that the characteristic struc- 
ture showing peak-dips in the noise spectrum may be 
changed after the electron-electron interaction is taken 
into account. More interesting phenomena in the noise 
spectrum should be expected in the interacting systems 
at low temperature. A closed formulation for the exact 
master equation of the reduced density matrix and the 
exact calculation of noise spectrum with the considera- 
tion of Coulomb interaction is not obvious. However, the 
present formalism is easy to be extended for including 
the Coulomb interaction with respect to the saddle point 
approximation^ or loop expansion;^ where the Coulomb 
interaction can be treated self-consistently in general- 
izing integrodifferential equations of motion Eq. Q, or 
by means of Hierarchical equation of motion approach. 
The work along this line is in progress. 
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Appendix: Numerical method 



As mentioned in Sec. Ill C[ the key to the practical cal- 
culation of the noise spectrum as well as the reduced den- 
sity matrix Eq. ^ and current in Eq. ^ in the present 
formulism relies on how to solve the integrodifferential 
equations of Eq. Q , which contain all the non-Markovian 
memory effects of the central system interacting with its 
environment. In general, it is not possible to analytically 
solve Eq. ^ but the numerical solution of such inte- 
grodifferential equation is also very difficult due to the 
non-Markovian memory kernels. Here, in terms of a pa- 
rameterized scheme^i^ we develop a numerical method 
in terms of a closed set of coupled differential equations 
of motion for solving u{t) and v{t) to overcome the dif- 
ficulties in the direct numerical calculation of the inte- 
grodifferential equations. 

For the sake of generality, we start with the energy 
dependence spectral density as a Lorentzian-type form 
centered at Ca^^^^^ 



r„(a;) 



(c^-e„)2 + W^2' 



(14) 



where Fq, describes the coupling strength and Wa is the 
line width of the source (drain) reservoir with a = L[R). 
Obviously the wide band limit (flat band), Ta{oj) = Fq,, 
is achieved by simply letting Wa oo. In terms of 
Eq. ([S]) , the non-local time correlation functions can be 
parameterized as22i2^ 



(15a) 



M 



M 



(15b) 



m=0 



m=0 



with gam{t ^ ''') = ^Qme~'*^°™'*~^^- The first term in 
da {t — t) with TO = arises from the pole of the spectral 
density function, with 

TaWal2 . 

The other terms with to > (M — > oo in principle) arise 
from the Matsubara poles, where the relevant parameters 
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are explicitly given as 



7q 



Pa 

(2m - 1)tt 



Pa 



(17a) 
(17b) 



The above paramterization is based on the Matsubara 
frequencies decomposition of Fermi distribution. An al- 
ternative efficient parameter scheme proposed in Ref.— 
is Fade spectrum decomposition for its convergence sig- 
nificantly faster than other schemes at all tempera- 
tures. The resulting formalisms of Eq. (jl5p are un- 
changed, but the coefficients for m > in g{t, r) are 
modified accordingly.'*^ Here, following the procedure in 
Refs. , 22.45. with the above parameterized scheme, we will 
present a numerical method in terms of a closed set of 
coupled differential equations of motion, instead of the 
integrodifferential equation of Eq. (|4ap , to solve the time- 
dependent transport current and the noise spectrum. 

Introducing a new function g'^{t) = J^* dr ga{t, t)u{t) 
which satisfies the following equation 

fifaW=9a(0)«(i)-7aoSr^(i), (18) 
with QaiO) — WaTa/2, then Eq. pa|) is reduced to 

uit)^-ie{t)u{t)-Y,9a{t)- (19) 

a 

In other words, the integrodifferential equation of 
Eq. (j4al) is transformed into a coupled pure differential 
equations of Eq. and Eq. ([TO]) , from which it is rather 
easy to obtain the time-dependent solution of u(t) and 
Qait)- On the other hand, the formal solution of Eq. (|4b| 
is 



(20) 



from which we have 



V{t) =v(i,t) = ^ dTi dT2u{t,Ti)ga{Ti,T2)u\t,T2). 

(21) 

Its time derivation is 

vit) = -i[e{t),v{t)] - J2 [^aW + H.c] (22) 



with Ia{t) — ICa{t) + Qa{t), where 



JCa{t) 



dr gait,T)v{T,t), 



(23a) 



Qa{t) = - / dTg^{t,T)u\T,t) ^ Qam{t). (23b) 



The second identity in Eq. (I23bp is based on the param- 
eterization scheme of gQ,(t,T) given by Eq. (|T5b)) . which 
leads to 

Qarn{t)=- f dr g^rn{t, t)u\t ,t) , (24) 

Furthermore, by introducing new functions 

(t) = dr dTiga{t,T)u{T,Ti)g {ri,t), 

Jto Jto 

(25a) 

'Daa'{t)= dr dTiga{t,T)v{T,Ti)ga'{Ti,t), (25b) 
Jto Jto 

we can find the following coupled differential equations: 

ICa{t) = ICo,it)[ie{t) - -f^o] + ga{0)v{t) 

+ [Caa'rnit) ~ V^^, [t)] , (26a) 

a'm 

Q.am{t) = Qarn{t)[ie{t) - -fam] - Varn + X!^a'am(0: 

(26b) 

Caa'mit) = -ffa (0) Q1,„ (t) - baO (0 + 7^'™ W]Caa'm (0 , 

(26c) 

Vc^a' {t) - ga (0)<, (t) + ICa {t)ga' (0) 
-[laQ{t)+ll,o{t)]'Daa'{t). 



(26d) 



Solving numerically the coupled pure differential equa- 
tion of Eq. (PB)) together with Eq. (1^ , we thus have the 
solution of v{t) and Ia{t) without directly calculating 
the multi-integrals in Eq. ([20)1 . Consequently, all the 
time-dependent coefficients in the transient current of 
Eq. ([6]) , the master equation of Eq. ([2]) and the n-resolved 
master equation of Eq. ([7]), as well as the frequency- 
dependent noise spectrum of Eq. can be easily ob- 
tained through the relations Ha{t) = g'^{t)[u{t)]~^ and 
\a{t) — 1-a{t) — fia{t)v{t). This numerical method 
largely simplifies the numerical difficulties in solving the 
integrodifferential equations of Eq. Q which involve very 
complicated non-Markovian memory kernels. 

In the wideband limit Wa oo, the numerical 
calculation can be largely simplified. Explicitly, the 
temperature-independent memory kernel of Eq. (I15ap is 
reduced to a delta function given by git.r) = ^5{t — 
t). For the temperature-dependent memory kernel in 
Eq. (jl5bp . the first term (m = 0) is also reduced to a 
delta function but the other terms (m > 1) are appar- 
ently changed not too much: 



ga{t - t) 



2(l + e-«/3=M/„) 

M 



m— 1 



(27) 



Thus the non-Markovian effect in the wideband limit is 
determined by the temperature together with the bias 
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voltage and the system-reservoir coupling. If we take 
further a high temperature limit /3q, — ^ 0, the summa- 
tion term in Eq. (f27| will also be reduced to a delta 
function of t — r. Then no any memory effect remains, 
and a true Markov limit is reached at high tempera- 
ture limit. In the practical numerical calculation, we 
keep the temperature-dependent memory kernel with the 
expression of Eq. (I15bp and the WBL is taken by set- 
ting Wa > 100r,Hii^ The equations of the m-related 
quantities, such as v{t), Qamit), and Caa'mit) given by 
Eq. dHI), Eg. ([26b)) and Eq. respectively, are thus 



not changed. The other equations can be simplified as 
follows. The equation of Eq. (IT51) is recast to 

r 

u{t) = ~ie{t)u{t) - -u{t), (28) 

and Eq. (IT5| is not needed. Also, we can simply solve 
Eq. ((23a|) and Eq. (I25bp with the results JCa = Tav{t)/2 
and Vaa'{t) = TaV{t)ra'/4: but Eq. (PBa)) and Eq. ([2M)) 
are also no longer used. 
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